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In the present work we analyze the second order transition line that connect the tricritical point 
and the quantum critical ending point on the temperature-magnetic-field plane in UGe 2 . For the 
microscopic modeling we employ the Anderson lattice model recently shown to provide a fairly 
complete description of the full magnetic phase diagram of UGe 2 including all the criticalities. The 
shape of the so-called tricritical wings, i.e. surfaces of the first-order transitions, previously reported 
by us to quantitatively agree with the experimental data, is investigated here with respect to the 
change of the total filling and the Lande factor for / electrons which can differ from the free electron 
value. The analysis of the total filling dependence demonstrates sensitivity of our prediction when 
the respective positions of the critical ending point at the metamagnetic transition and tricritical 
point are mismatched as compared to the experiment. 
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I. MOTIVATION AND OVERVIEW 

Quantum critical phenomena have captured general at¬ 
tention due to their unique singular properties observed 
at low temperature (T ->• 0) and near the quantum criti¬ 
cal point (QCP) which is frequently accompanied by the 
unconventional superconductivity (SC) [T]. From this 
perspective, /-electron compound UGe 2 is a system with 
phase diagram comprising coexistence of spin-triplet SC 
and ferromagnetism (FM) |2HS|> as well as an abundance 
of critical points (CP), either of quantum and classical 
nature |7]. Experimental studies among others have re¬ 
vealed existence of the two characteristic classical CPs in 
the absence of the field (cf. Fig. [^: (i) the critical end¬ 
ing point (CEP) at 7K at the metamagnetic transition 
separating strong (FM2) and weak magnetization (FMl) 
regions iBHin], and (ii) the tricritical point (TCP) at the 
FM to paramagnetic (PM) phase transition located at 
T = 24 K. Additionally, with the applied magnetic field 
the second order transition line starting from TCP can 
be followed to T = 0 where it is expected to terminate in 
a quantum critical ending point (QCEP) OHO]. In effect 
the magnetic phase boundaries in UGe 2 reflects so-called 
tricritical wing shape. 

Such complex magnetic phase diagram with all the 
above criticalities, both classical and quantum, is partic¬ 
ularly challenging in terms of theoretical modeling. One 
of the first approaches, based on the single-band model 
describing tricritical wings, was the work by Belitz, et al. 
HU. However, the microscopic description of the mag¬ 
netic phase diagram with all the CPs including also CEP 
at the metamagnetic transition, as observed in UGe 2 , has 
been missing until our recent works Hsmu. 

Our analysis is based on the (two-orbital) Anderson 
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Figure 1. Schematic magnetic phase diagram of UGe 2 on 
pressure-temperature plane drawn on the basis of the exper¬ 
imental results [5]. 


lattice model (ALM) [131 HI], often referred to as the peri¬ 
odic Anderson model. Findings for UGe 2 , both from first 
principle calculations and experiments are: the quasi- 
two-dimensional character of the Fermi surface HU, a 
uniaxial anisotropy for magnetization |16| . U-U inter¬ 
atomic distance above the so-called Hill limit (U > f^nd the 
paramagnetic moment per U atom different from that ex¬ 
pected for either or atomic configurations Emu. 
We show in the following that all of these findings can be 
coherently explained within our two-orbital model start¬ 
ing with originally localized /-states and subsequently 
being strongly hybridized with the conduction (c) band 
states on a two dimensional lattice and with the applied 
magnetic field accounted for by the Zeeman term only. 

Ferromagnetic order in our model arises from effect 
of competing hybridization and the /-/ interatomic 
Coulomb repulsion. The emergence of two distinct ferro¬ 
magnetic phases is in our model driven by the changing 
topology of the Fermi surface HM which in turn is 
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induced by a relative motion of hybridized and spin split 
subbands with the increasing f-c hybridization. The re¬ 
sults obtained from such picture m qualitatively agree 
with the majority of UGe 2 magnetic and electronic prop¬ 
erties, as seen in neutron scattering m, de Haas-van 
Alphen oscillations p2ll2^ . and magnetization measure¬ 
ments [S]. Also, a semi-metallic character of the weak 
FMl phase is supported by the band-structure calcula¬ 
tions m- A similar idea concerning the emergence of 
two distinct FM phases in UGe 2 was also obtained earlier 
within the phenomenological picture based on the Stoner 
theory incorporating a two-peak structure of the density 
of states in a single band [^. In brief, our microscopic 
model extended to the case of T > 0 [H] describes well 
emergence of all GPs on the magnetic phase diagram of 
UGe 2 |5HTD] in the semiquantitative manner m- Here 
we compare in detail our results with the experimental 
data, namely predicted second order transition line join¬ 
ing the TGP with the corresponding QGEP. In particular, 
we determine the influence of the following factors: (i) 
the total band Ailing n, and (ii) the value of the Lande 
factor gf for f states, on the position of this second-order 
line. The influence of factor (i) has the following impor¬ 
tance. For exemplary Ailing n = 1.6 we have shown I Ml 
that the relative position of the TGP and GEP (cf. FigQ 
is the same as that seen in the experiments UdD]. The 
important question is whether such a mutual alignment 
of those two critical points is necessary to achieve a good 
fit and to what extent the proper curvature of the line 
joining TGP and QGEP is robust with respect to the se¬ 
lected band Ailing. The discussion of the dependence on 
(ii) has its justification in the not fully resolved nature 
of magnetism in heavy-fermion systems in general and 
UGe 2 in particular. Although it is assumed and widely 
accepted to be fully itinerant [5], there is evidence for a 
partially localized contribution |24l EH] . In such case, the 
influence of the orbital effects and their coupling to the 
spin should have an influence on gj value. 


II. MODEL AND APPROACH 

We begin with the orbitally nondegenerate Anderson- 
latice model (ALM) on square lattice and with applied 
magnetic field accounted for via the Zeeman splitting 
(i.e., with the effective field is h = ^ggofisH), so that 
the starting Hamiltonian is 

"Ho = E fi-^ 

ij,t7 i,a 

+ + ( 1 ) 
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where the onsite hybridization is of magnitude P < 0 and 
the Lande factor for / electrons is gf (the free electron 
value is 5 = 2). The model describes a two-orbital system 


with the conduction (c) band arising from the nearest 
(t) and the second-nearest (t') neighbor hoppings, and 
the strong /-/ Goulomb interaction is of magnitude U. 
If it is not stated otherwise, we set t' = 0.25|t|, U = 5|t|, 

£/ = 5 / = S = 2, and n E Z^(hf^ + n{^) = 1.6. 

We also add to the Hamiltonian 0 the usual term 
with the chemical potential /i, i.e., 

= + ( 2 ) 

i,cr 

The model is solved here by means of statistically con¬ 
sistent Gutzwiller approximation (SGA) |27II29] . The 
method was successfully applied to the number of prob¬ 
lems [sniEi]. It is characterized with the physical trans¬ 
parency and flexibility, that it could be also incorporated 
into other methods such as EDABI |32l 133] . 

We introduce the Gutzwiller projection acting onto un¬ 
correlated wave function |i/)o) in the following manner 

\^g) = UPG;M, (3) 

i 

where |'0g) is the wave function of the correlated ground 
state. In effect, we map many-particle correlated Hamil¬ 
tonian (Q onto an effective single-particle Hamiltonian 
T~Lsga acting on uncorrelated wave function |V'o)) that, 
after taking the space Fourier transform, is as follows 


iisGA = 


An/)-A^( E <^nL- 
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( 4 ) 


where e (Ck^j/ko-)) A is the number of the sys¬ 
tem sites, Qa- is the hybridization narrowing factor, and 
(f = )o m- Necessary constraints for the / elec¬ 

tron number and their magnetic moment |29) are incorpo¬ 
rated by means of the Lagrange multipliers and A^, 
respectively. Hamiltonian 0 can be straightforwardly 
diagonalized with the resulting four eigenvalues {E^^} la¬ 
beled with the spin (cr = ±1) and hybridized-band (6 = ±1) 
indices. For T > 0 we construct a generalized Landau 
grand-potential functional according to 
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which is next minimized with respect to the set 
of following parameters assembled to a vector, A e 
{d, n/,TO/, A^, A^}. Additionally, we adjust self- 
consistently the chemical potential from the con¬ 
dition of fixing the total number of electrons, 
n = l/AX!k&(T/(^kcr)i where f{E) is the Fermi-Dirac 
function. Finally the ground state energy is defined by, 


Eg = 


( 6 ) 
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where means that the value of T is taken at minimum 
for the parameters A. 

III. RESULTS AND DISCUSSION 

We start our analysis with the discussion on the proper 
assignment of the physical units to the microscopic pa¬ 
rameters provided so far in dimensionless units (i.e. 
scaled by |t|) to make the quantitative comparison with 
the observed UGe 2 characteristics. To do so, we have 
adjusted them [14] by matching the relative positions of 
the two classical CPs: TCP and CEP at the ferromag¬ 
netic transition, as well as attributing the experimentally 
measured critical temperatures. Matching the results in 
physical units by fixing the position of two critical points 
we would call a strong fitting, whereas by fixing the po¬ 
sition of just a single of them a weak fitting. 

In our previous works [min] we have found that for 
the total filling n = 1.6 we could coherently and quantita¬ 
tively describe the UGe 2 phase diagram. Although, there 
is no direct experimental evidence in UGe 2 for choos¬ 
ing this particular filling, we have first matched for cho¬ 
sen n TCP and CEP temperatures according to the ex¬ 
periment na - strong fitting condition, and second, we 
have verified our prediction obtaining agreement of the 
second-order transition line joining TCP with QCEP m 
with that measured. Additionally, the comparison has 
provided among others the estimate of the QCEP ap¬ 
pearance about 30 T, i.e., higher than that suggested in 
Ref. [TU| which is 18 T. 

A natural question arises if this test is sensitive to the 
choice of n. We show in Fig. that the second-order 
transition line joining TCP and QCEP determined for a 
slightly different total filling than n = 1.6 deviates sig¬ 
nificantly from the trend of the experimental data m- 
Hence indeed, the comparison is very sensitive to the 
choice of n. Thus, together with equally sensitive adjust¬ 
ment of Ttcp and Tqep with respect to the choice of 
n m, it is unlikely that our excellent agreement for the 
single value of parameter n is fortuitous. 

For the sake of completeness and reference to other re¬ 
lated works we include in Fig. the (dashed) curve pre¬ 
dicted by the one of the most successful approaches de¬ 
scribing the general tricritical behavior in itinerant mag¬ 
netic systems mi HD. In that procedure, the necessary 
inputs are the positions of the two CPs, namely TCP 
and QCEP, leading in fact to the strong fitting, but with 
different pair of CPs. However, such fitting can be as¬ 
sociated with an error as the position of the QCEP, in 
contrast to TCP and CEP, is not experimentally deter¬ 
mined but only extrapolated to 18 T, following Ref. |10| . 
Note that this condition in our modeling is satisfied for 
the total filling around n 1.55 (cf. Fig. if we take 
the extrapolated value of the critical field. In this case, 
the comparison with results for UGe 2 [10] is worse then 
e.g. for n = 1.6 and the temperature of the CEP is much 
lower than that determined in experiments [^. 



Figure 2. The second-order transition line joining TCP and 
QCEP for selected values of total filling, n. Solid, red curve is 
replotted from the results reported in Ref. mi. Experimental 
points are extracted from pictures in Ref. m- Grey, dashed 
curve predicted by the theory formulated by Belitz, et al. m 
is shown for comparison (see main text). 


It is worth mentioning that the model employed by 
us belongs to the class discussed earlier by Kirkpatrick 
and Belitz m to reflect the generic tricriticality in the 
case of metallic magnets. Namely, systems in which the 
conduction electrons are not a source of the magnetism 
themselves, but are couple to the magnetically ordered 
localized electrons in a second band. The origin of the 
first-order transition at low temperature described within 
the mean-field theory developed in the Refs. mm is 
based on the effect of the soft fermionic modes coupled to 
the magnetization fluctuations, and thus differs from our 
approach. Here the mechanism for ferromagnetism is due 
to the coupling of the conduction electrons with localized 
/ states by hybridizing with them and competing with 
the /-/ Coulomb interaction. This competition in the 
Stoner-like manner induces phase transitions associated 
with the abrupt changes of the Fermi surface topology. 

The simplest verification of our analysis can be carried 
out by means of chemical alloying, i.e., by changing the 
electron concentration in the system. However, the lack 
of known isostructural compounds to UGe 2 may be an 
apparent obstacle for such test. Though, the determina¬ 
tion of the tail of the 2nd order line joining TCP and 
QCEP for the field larger than 16 T should provide an 
insight on the issue whether our model correctly predicts 
the appearance of QCEP around 30 T m- 

If our model is to be used to understand the magnetism 
of other ferromagnetic superconductors: URhCe [M] and 
UCoGe [3S], it would provide a perfect testing ground 
of our model as those compounds have been frequently 
studied by means of the chemical substitution [IJI3H6]. 

Finally, we provide a brief analysis of the impact of 
the Lande factor value for / electrons, gj, i.e., in the 
situation when the 2 component of the total spin of the 
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Figure 3. The second-order transition line joining TCP and 
QCEP for selected values of Lande factor for / electrons, gf. 
Solid curve for gf = 2 is replotted from the results reported in 
Ref. |14| . Experimental points are extracted from Ref. |10| . 


system does not commute with Hamiltonian. In Fig. 
we present the curves for three different values of gf. 
The curve for gf = 2, is plotted as the reference curve 
and is based on the results of Ref. M- Value of gf 
is not known for UGe 2 and generally, for complex com¬ 
pounds is has a tensor character which depends on the 
magnitude of the spin-orbit coupling. For that reason 
we restrict our discussion to the comparison when g/ is 
equal to the free electron value gf = 2, and subsequently 
when is lower and higher (cf. Fig® [36]. Specifically, the 
lower value of Lande factor = 6/7 is motivated by that 
for the Ce-based compounds, where it can be derived for 
the spin S' = 1/2 and angular momentum L = 3, oriented 
antiparallel and where, strictly speaking, our model is 
also generally valid, as long as we do not account for the 
orbital degeneracy of / states of the uranium-based ma¬ 
terials. As presented in Fig. it seems that any value of 
g factor for / states which deviates considerably from the 
free-electron value, provides much worse agreement with 
experimental data m- In conclusion, due to predomi¬ 
nantly itinerant nature of / electrons in UGe 2 jH], it is 
very likely that any crystal-field derived multiplet struc¬ 
ture is washed out and hence the value <?/ - 2 should be 
regarded as realistic value. Nevertheless, the problem of 
double localized-itinerant nature of /-electrons |24l [26] 
may arise as the system evolves with the increasing tem¬ 
perature, in comparison to the pressure evolution at low 
temperature studied in detail here. 

IV. OUTLOOK 

In the present work we have employed the Anderson 
lattice model HllTlj to provide a fairly complete descrip¬ 


tion of the magnetic phase diagram (p-T-h profiles) of 
UGe 2 including all the criticalities for this compound. In 
particular, we study the effect of the choice of the total 
filling on the quality of the fit, based on our model, to 
the experimental data m concerning the second-order 
transition line joining the critical points TGP and QGEP. 
We have found that our prediction is very sensitive to the 
change of n, which leads also to a mismatch of critical 
temperatures of TGP and GEP at at the metamagnetic 
transition as compared to the experiment. We infer from 
this result that our excellent agreement for the single 
value of n is unlikely fortuitous. We have also analyzed 
the effect of the Lande factor gf value for / electrons. 
In this case, any sizable deviation from the free electron 
value gf = 2 shifts the theoretical curves away from the 
experimental points. Thus treating / electrons as truly 
itinerant electrons in UGe 2 seems to be fully justified. 

Our final remark addresses the problem of the spin- 
triplet superconductivity (SG) origin occurring in UGe 2 
in [5] . We have predicted in our previous work [Tl] the 
appearance of QGP in the vicinity of the SG dome. It 
have been proposed that that GEP (cf. Fig. Q at the 
metamagnetic phase boundary can be followed down to 
the T = 0 by changing both the electron concentration 
and the hybridization magnitude \V\ (cf. Fig. |^. The 
proposed quantum critical point is of Lifshitz type as it 
separates states with two distinct Fermi-surface topolo¬ 
gies. Quantum critical fluctuations or the residual /-/ 
Hund’s rule interaction (neglected here) can become the 
possible source of the spin-triplet superconductivity |27|- 
l4T] . A detailed and quantitative discussion of the pairing 
requires a separate analysis. 
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Figure 4. Evolution of the critical temperature of CEP to¬ 
wards Lifshitz QCP driven by changing total filling, n and 
hybridization, V reproduced after mi. 
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